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ABSTRACT 

It is well established that timing and spectral properties of Galactic Black Hole 
(BH) X-ray binaries (XRB) are strongly correlated. In particular, it has been 
shown that low frequency Quasi-Periodic Oscillation (QPO) vi ow - photon index 
T correlation curves have a specific pattern. In a number of the sources studied 
the shape of the index- low frequency QPO correlations are self-similar with a 
position offset in the — P plane determined by a BH mass M B h • Specifically, 
Titarchuk & Fiorito (2004) gave strong theoretical and observational arguments 
that the QPO frequency values in this ui ow — P correlation should be inversely 
proportional to M B h- A simple translation of the correlation for a given source 
along frequency axis leads to the observed correlation for another source. As a 
result of this translation one can obtain a scaling factor which is simply a BH 
mass ratio for these particular sources. This property of the correlations offers 
a fundamentally new method for BH mass determination in XRBs. Here we 
use the observed QPO-index correlations observed in three BH sources: GRO 
J1655-40, GRS 1915+105 and Cyg X-l. The BH mass of (6.3 ±0.5) M 0 in GRO 
J1655-40 is obtained using optical observations. RXTE observations during the 
recent 2005 outburst yielded sufficient data to establish the correlation pattern 
during both rise and decay of the event. We use GRO J1655-40 as a standard 
reference source to measure the BH mass in Cyg X-l. We also revisit the GRS 
1915±105 data as a further test of our scaling method. We obtain the BH mass 
in Cyg X-l in the range 7. 6-9. 9. 

Subject headings: accretion, accretion disks — black hole physics — starsdndividual 
(Cyg X-l), individual (GRO J1655-40), individual (GRS 1915-1-105) :radiation 
mechanisms: nonthermal — physical data and processes 
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1. Introduction 

In Galactic Black Hole (GBH) sources transitions between spectral states are accompa- 
nied by a characteristic evolution of timing properties (see McClintock & Remillard 2006, 
for the review of BH spectral states). In particular, Low Frequency Quasi-Periodic Oscilla- 
tions (LF QPO) observed during low-hard and intermediate states are very closely correlated 
with the photon index of the power law spectral component (e.g. Vignarca et al. 2003). 
Recently, the same types of correlation have been firmly established for Cyg X-l [Shaposh- 
nikov & Titarchuk (2006), hereafter ST06] and GRO J1655-40 [Shaposhnikov et al. (2007), 
hereafter S07] . In this Letter we use these correlations to constrain the mass of BH in Cygnus 
X-l based on the theoretically motivated dependence of QPO-index pattern on the mass of 
the central BH. 

Some progress has been made in understanding how variability scales with BH mass 
between GBHs and AGNs (see Done, & Gierlinski, 2003, and references therein). However, 
in most works authors concentrate on the observational aspects giving less attention to the 
physical reasoning behind the scaling method. The ubiquitous nature of the correlation 
of index with QPO frequency suggests that the underlying physical process or condition 
which gives rise to the LF QPO is tied to the corona; and, furthermore, that this process 
varies in a well defined manner as the source progresses from one spectral state to another. 
Moreover, the fact that the same correlations are seen in so many sources, which vary 
widely in both luminosity (presumably with mass accretion rate) and state, suggests that 
the physical conditions controlling the index and the low frequency QPOs are characteristics 
of these sources, and that by virtue of the low-high-frequency correlations [see Psaltis, Belloni 
& van der Klis (1999) and Belloni, Psaltis & van der Klis 2002), they may be a universal 
property of all accreting compact systems. The observational evidence presented above 
motivated to develop a detailed physical model of the corona surrounding a BH which directly 
predicts the behavior of the spectral index with fundamental properties of the corona [see 
Titarchuk, Lapidus & Muslimov (1998), hereafter TLM98; Laurent & Titarchuk (1999); 
Titarchuk & Osherovich (1999) and Titarchuk & Fiorito (2004), hereafter TF04). This 
model incorporates fundamental principles of fluid mechanics, radiative transfer theory and 
oscillatory processes. It identifies the origin of the QPO as a fundamental property of a 
compact coronal region near the BH and shows how the photon index of this corona changes 
as a function of mass accretion rate. It is important to emphasize that the variable QPO 
frequency scaling as 1/M is a generic feature of the TLM98 model. In fact, the QPO 
frequency scales as a ratio of plasma (sound) velocity and the size of the oscillating region, 
which is measured in the Schwarzchild radius units. 

TF04 illustrate how the mass of one galactic black hole (BH) source can be used to 
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determine that of another using the RXTE data from GRS 1915+105 and XTE J1550-564, 
which exhibit remarkably similar QPO frequency— index correlation curves. The scaling 
factor between the LF QPOs of two X-ray BH sources with similar power-law (PL) photon 
indices can be used to determine the ratio of their BH masses. The mass determination using 
the QPO frequency— index correlation is consistent with X-ray spectroscopic and dynamical 
mass determinations for these sources (see references in Shrader & Titarchuk 2003, hereafter 
ST03). The scaling was successfully applied for ultra- luminous sources using QPO frequen- 
cies (Fiorito & Titarchuk 2004; Strohmayer et al. 2006). Fiorito & Titarchuk (2004) and 
Dewangan, Titarchuk & Griffiths (2006) applied this new method of BH mass determination 
to estimate the black hole mass of the ultraluminous X-ray sources M82 X-l. Using scaling 
arguments and the correlation derived from the consideration of Galactic black holes, they 
conclude that M82 X-l is an intermediate black hole with a mass in the range of (250-1000) 
M 0 . Recently this technique led Strohmayer et al. (2006) to estimation of BH mass of 
~ 1500 Mq in ULX source NGC 5408 X-l. 

In this Letter we demonstrate how the TF04 scaling method can be applied to the 
determination of the Cyg X-l BH mass M Cyg x- 1 - We obtain Mc yg x - i when we slide the 
index-QPO correlation of GRO J1655-40 vs that of Cyg X-l along frequency axis. We 
test this scaling method using BH masses of GRS 1915+105 and GRO J1655-40 previously 
inferred from optical and IR observations. 


2. Sources and Observations 
2.1. GRO J1655-40 

The X-ray binary GRO J1655-40 (GRO J1655 hereafter) is a well-known example of a 
BH X-ray transient that underwent several major outbursts within the last 12 years. It was 
discovered by the BATSE instrument onboard the Compton Gamma- Ray Observatory in mid 
1994 (Zhang et al. 1994). The secondary star being relatively bright, the binary parameters 
are exceptionally well determined among LMXBs. The most recent optical photometry of 
Green, Bailyn & Orosz (2001) led to a BH mass estimate of (6.3 ± 0.5 )M©. Hjellming & 
Rupen (1995) inferred a distance to GRO J1655 of (3.2 ± 0.2) kpc from radio jet analysis. 
During RXTE era GRO J1655 had two major outburst in 1996/97 and in 2005. However, 
during the first event pointed RXTE observations started when the source was already in the 
high-soft state and the state transition was not covered. The 2005 outburst was recognized 
at a very early stage and observers at different wavelengths were relatively prepared and 
successful at responding relatively quickly (see S07, for extended report on multiwavelength 
campaign during 2005 outburst). Almost daily RXTE coverage allowed to track down a 
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timing and spectral evolution of the source in detail. We distinguish the data during the 
outburst rise and decay because the correlation patterns are slightly offset presumably due 
to different plasma temperatures in the corona (see below for explanation). The data during 
the outburst rise were collected between MJD 53420 and 53440, while the decay data were 
collected between MJD 53627 and 53640. We refer a reader to S07 for the details of spectral 
and timing data analysis of the outburst rise data. The decay data were analyzed in the 
same manner. 


2.2. Cyg X-l 

Cyg X-l is one of the brightest high-energy sources in the sky, with an average 1- 
200 keV energy flux of ~ 3 x 10 -8 ergs cm -2 s -1 . Its optical companion is an 09.7 lab 
supergiant HDE 226868. Estimates of the mass of the X-ray star, 5^M 0 <15 [e.g., Herrero 
et al. (1995)] strongly suggest the presence of a black hole. Observed spectral and temporal 
X-ray characteristics are extensively studied by ST06, based on the large amount of data 
collected in the RXTE archive (see §2 in ST06 for the data description). They presented 
the analysis which includes ~ 2.2 Ms of RXTE archival data from Cygnus X-l. The main 
findings of the ST06 study are following: The photon index T steadily increases from 1.5 
in the low/hard state to values exceeding 2.1 in the high-soft state. The low frequency vl 
is detected throughout the low-hard and intermediate state, while it disappears when the 
source undergoes transition to the soft (thermal dominated) state. Like in other BH sources, 
there is an indication of saturation of the index in the high-soft state for Cyg X-l (see Fig. 
1). This saturation effect, which is presumably due to photon trapping in the converging 
flow, can be considered to be a BH signature. 


2.3. GRS 1915+105 

GRS 1915+105 (GRS 1915 hereafter) is another prominent Galactic BH binary. The 
source is very bright and exhibits very diverse temporal properties. In attempt to classify 
the variability in GRS 1915 Belloni et al. (2000) identified at least 12 different classes. 
Morgan, Remillard & Greiner (1997) reported three different types of QPO: a QPO with 
constant centroid frequency of 67 Hz, a dynamic low-frequency (10~ 3 — 10 Hz) QPO with 
a large variety of amplitudes and widths, and high- amplitude “sputters” at frequencies of 
10 -3 — 10 -1 Hz. Authors also mention that it is difficult to describe the source properties 
in term of conventional BH states. A BH mass estimation for GRS 1915+105 of 13.3 ± 4 is 
available from IR observations (Greiner et al. (2001)] and X-ray observations [Borozdin et 
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al. (1999) and Shrader & Titarch.uk (2003)]. 

Despite the complicated nature of the GRS 1915 spectral and timing properties we were 
able to construct the QPO-photon index correlations. We use the data from six RXTE ob- 
servations, namely: 10408-01-27-00, 20402-01-50-00, 20402-01-50-01, 20402-01-51-00, 20402- 
01-52-000, 40702-01-02-00. All of these observations, except the last one, axe taken during 
the hard state. During this state a power spectral density (PDS) of the source is very close 
to PDS of classical hard state (% Belloni variability class), i.e. the band-limited noise is 
accompanied by one or two low frequency QPOs. Data from these observation yield the 
lower part of QPO-index correlation curve (see Fig. 2). During observation 40702-01-02-00 
the exhibited irregular variability with flaring episodes (v class). This observation yields the 
upper part of the correlations with a saturation. For each observation we calculated energy 
and Fourier spectra for consecutive intervals of 128 seconds. Following the approach of Vi- 
gnarca et al. (2003) we fit the data with a sum of multicolor disk and power law models. 
We fit PDS continuum with broken power law shape to model band-limited noise and we 
modeled QPOs with Lorentzians. 


3. QPO-spectral index correlation model 

In most BH sources, the correlation of the spectral index— QPO frequency follows a 
remarkably similar pattern. Namely, at lower QPO frequencies the dependence starts at the 
index of ~ 1.5 (e.g. low-hard state) and follows linear correlation with positive slope until 
some QPO frequency value v br where it levels off. This saturation of index versus QPO 
frequency approximately corresponds to the intermediate state. Then, upon the QPO turn 
off, the source enters the soft state. 

According to our scenario (see TLM98 and TF04) the LF QPO features are observational 
appearances of a Compton cloud (CC) undergoing normal mode volume oscillations. The 
radius of the CC outer boundary is a radius of an adjustment of disk Keplerian rotation to 
sub-Keplerian CC rotation should be determined in Schwarzchild radius (Rs) units, namely 
Rout = x out Rs ■ Thus the CC size, as (R out — ) should be proportional to m = M/M Q 

(Ri n is also scaled with Rs) and QPO frequency (as a ratio of the velocity and size) should 
be reciprocal to the characteristic size of the system and consequently to m. The spectral 
index, on the other hand, is a measure of CC upscattering efficiency. 

Note that direct connection of variability in the low-hard and intermediate states with 
power law spectral component has been identified and well established in many studies (see 
van der Klis 1995, and references therein). In the hard state the mass accretion rate is low, 
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the size of the corona is large and the index is determined by thermal Comptonization in an 
optically thin geometrically thick configuration. When accretion rate starts to increase a cold 
disk becomes stronger and CC cooling by disk photons becomes more efficient. Ultimately 
the CC plasma temperature goes down. In the high-soft state the CC temperature is of 
the order of the disk temperature (about 1 keV). As the mass accretion rate increases the 
adjustment radius decreases and consequently the CC cloud contracts (see TLM98, T099). 
Thus QPO frequency increases with the CC contraction, while decreasing CC temperature 
leads to the lower efficiency of Comptonization, and thus to higher photon indices. 

The presented picture explains a positive correlation between spectral index and QPO 
frequency. The index saturation for higher QPO values can be accounted for by the bulk 
motion inflow (dynamical) Comptonization which becomes dominant for the final stage of 
state transition when the corona is cooled down by the disk photons. The saturation level 
of index can vary for the same source for different transition episodes. This index saturation 
value is determined by plasma temperature of the converging flow (see Laurent & Titarchuk 
(1999) during a transition. Details of this physical scenario theory can be found in Titarchuk, 
& Fiorito (2004) (and references therein). 

We must stress that the main property of index-QPO relation which we use for scal- 
ing method does not depend on a specific physical mechanism of photon and power density 
spectral formation. Namely, we only use that QPO frequency is inversely proportional to 
the size of the Compton cloud [i.e. corona, converging flow (ADAF) region] which, in turn, 
scales with the Schwarzshild radius and thus it scales with a BH mass (see TLM98, TF04). 
Moreover, the other important advantages of the technique is its independence of the system 
geometry and a distance to a source. The downside of the method is that the mass of one 
BH exhibiting clear index-QPO pattern has to be well known. Fortunately, such a system is 
available, i.e. GRO J1655-40. 


4. Data analysis. Application of scaling method 

In principle, it is possible to construct index-QPO relation on the basis of the physical 
theory presented in TF04. However, the exact relation between optical depth r and Reynolds 
number 7 (which is inverse of a-parameter in the disk) requires the disk equation of state. 
The direct implementation of the theory has to include some assumption about the functional 
form of r( 7) because of uncertainties of the disk equation of state and the plasma viscosity 
in the disk. It is worth noting that any a— disk model does not specify an exact value of 
a— parameter [see e.g. Shakura & Sunyaev (1973); Narayan & Yi (1994); Chakrabarti & 
Titarchuk ( 1995)] . The fact that this scaling method follows from- -the first principles the 
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QPO-BH mass scaling allows us to avoid the full implementation of the theory to construct 
QPO-index relation for each particular source. All we need is a simple function which 
correctly represens the data and give us means to calculate a scaling coefficient between two 
correlation curves. The correlations between QPO frequency and photon index for three 
sources considered in this Letter are shown on Figures 1 and 2. It is clear that QPO-index 
correlation follows a simple pattern. Namely, the dependence is linear with particular slope 
for lower frequencies, while near some frequency (which we call v tr ) the function smoothly 
levels off to become constant. This behavior is well reproduced by the following analytical 
function: 

f(v) = A - D£ln{exp[(i/ ir - v)/D] + 1}. (1) 

The two asymptotes of this function are f(y ) = A + B(v — v tr ) for u < v tr and f(u) = A for 
v > u tr . Thus, A is a value of the index saturation level, B is a slope of the low frequency 
part of the data, u tr is the frequency where at which index-QPO dependence levels off . 
Parameter D controls how fast the transition occurs and it was fixed at 1.0 Hz throughout 
the fitting procedure presented in the next section. It is important to point out that the 
slope of B of the asymptote for v < v tr in Eq. (1) strongly depends on the BH mass value. 
Using the above model function we fit the observed index-QPO correlations and apply the 
scaling method for mass BH determination. 


4.1. Scaling from GRO J1655-40 correlation to GRS 1915+105 correlation 

For scaling with GRS 1915 we choose the data from the rise of the 2005 GRO J1655 
outburst because in both cases the photon index saturation level is approximately T ~ 2.3. 
We first fit the GRO J1655 points. We found the best fit parameter values of: Aji655 = 
2.28+0.02, .8/1655 = 0.133+0.008 Hz -1 and r+jiess = 6.64+0.48 Hz. The D was fixed at 1.0 
Hz. The best-fit curve is shown by a red solid line in Figure 2. The best fit values found for 
GRS 1915+105 are: A1915 = 2.25 i 0.01, R19X5 — 0.33 i 0.02 Hz - * and i^tr , 1915 = 2.23 i 0.07 
Hz. For a scaling factor we have s t 7i 655 _ +1 9i5 = ! R1915/R/1655 = 2.57 ± 0.17. We obtain a BH 
mass in GRS 1915 as 

Afi9i5 = Rl915 x Afj!655/Rji655 = (12.8 — 19.0)M©, 

using the scaling method and the value of BH mass in GRO J1655 of (6.3 ± 0.5 )M©. The 
inferred uncertainties (error bars) of Mi 915 is mostly affected by scattering and error bars 
of index-QPO points in the GRS 1915 correlation (see Fig. 2). The value of BH mass of 
(15.9+3.1 )M© is consistent with the previous BH mass estimate in GRS 1915 of (13.3±4)M 0 . 
This result of QPO-index correlation scaling demonstrates a potential power of the new 
method of a BH mass determination. 
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In Figure 3 we illustrate the idea of the scaling method using individual PDSs related 
to the same spectral state (index). The PDSs presented for GRO J1655 and GRS 1915 have 
the same flat-top noise levels and similar slopes after break but QPOs for GRO 1655-40 are 
shifted by approximately a factor of ~ 2.5. 


5. Scaling from GRO 1655 to Cyg X-l correlations. BH Mass in Cyg X-l. 

Now we proceed with the main result of this Letter, namely, with the mass determination 
of BH in Cyg X-l using QPO-index scaling. We apply exactly the same procedure (described 
in the previous section) for the RXTE mission-long data of Cyg X-l and the GRO J1655 
data from decay stage of its 2005 outburst. We perform the initial fit to Cyg X-l data to 
obtain: Ac yg x-i = 2.00 ± 0.01, Bc yg x- i = 0.065 ± 0.001 Hz -1 and Vtr,c yg x-i = 7.68 ± 0.15 
Hz. Fit for GRO J1655 decay data gives: Aji &55 = 2.05 ± 0.02 and Bji &55 = 0.047 ± 0.002 
Hz -1 and ^ r ,ji 655 = 11.4 ± 0.6 Hz. For the mass of BH in Cyg X-l we have: 

Mc yg x-i = (7.6 — 9.9 )Mq. 

Here we again use (6.3 ± 0.5) M© as a value of BH mass in GRO J1655. This inferred value 
of Cyg X-l BH mass is slightly lower than 10.1 M© given by Herrero et al. (1995) (with no 
error bars quoted) and it is within the mass range found by Gies & Bolton (1986). 

In Figure 4 we present the individual PDS of Cyg X-l and GRO J1655 related to the 
same spectral state (index). One can see the corresponding values of QPO frequencies are 
very close. However QPO in Cyg X-l is less coherent than that in GRO 1655 which is 
presumably a result of photon scattering in the strong ambient wind of Cyg X-l (see ST06 
for details). 

6. Conclusions 

Our results for the QPO-index correlation scaling from the GRO J1655 to GRS 1915 
and from the GRO J1655 to Cyg X-l data points demonstrates a potential power of the new 
method of a BH mass determination. 

We have tested this scaling method using the known (from optical and IR observations) 
BH masses in GRO J1655 and GRS 1915. We were able to reproduce and confirm the 
previous estimate of BH mass of GRS 1915 M igi5 . We obtain that M 1915 /M© = 15.9 ± 3.1. 

Our scaling method along with the observed index-QPO correlation give for the BH 
mass in Cyg X-l Mc yg x- i = (7.6 — 9.9 )M©, which is the most close constrains on the BH 
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mass for this prominent source obtained up to date in literature. Low Cyg X-l mass function 
is an obstacle for dynamical mass measurement from optical data, which also affected by 
the uncertainty in the source distance. Our alternative method does not suffer from these 
restrictions, which allows such a precise measurement of the BH mass. 

We acknowledge the referee’s constructive suggestions for the improvement of the paper 
presentation. 
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Fig. 1. — Determination of BH mass in Cyg X-l. Correlations for Cyg X-l comprises 
RXTE-mission-long data from ST06 (blue color). Data from GRO 1655-40 is for the 2005 
outburst decay. These sources show similar high frequency saturation levels for each of the 
correlations. Index-QPO correlation scaling gives us the BH mass in Cyg X-l in the range 
of (7.6 -9,9) M e . 
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Fig. 2. — Spectral Index-QPO frequency correlations data for GRS 1915+105 (blue points) 
and GRO J1655-40 during the 2005 outburst rise (red points). A clear saturation of index for 
high values of QPO is seen for both sources. The lower part of correlation is not available for 
GRS 1915+105. The data are fitted with with a model function to apply a scaling method. 
We obtain the BH mass in GRS 1915+105 in the range of (12.8 — 19.0) Mq using the scaling 
method and the well-known BH mass of GRO 1655-40 which is in the range of 6.3 ± 0.5 Mq. 
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Fig. 3.— Individual PDS for GRS 1915+105 (RXTE Observation ID 20402-01-50-00, 128 
sec. interval, see text for details) and GRO J1655-40 (RXTE Observation ID 90704-04- 
01-00) illustrating how low-frequency QPOs scale with the mass of the central black hole. 
The photon indices are 1.68 ± 0.02 and 1.67 ± 0.08 correspondingly. PDS have the same 
flat-top noise levels and similar slopes after break but QPOs for GRO 1655-40 are shifted by 
approximately a factor of 2.5. 
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Fig. 4. — Individual PDSs for Cyg X-l (RXTE Observation ID 50119-01-03-01) and GRO 
1655-40 (RXTE Observation ID 91702-01-79-00). Spectral index for both observations is 
about 1.92. QPO feature in Cyg X-l is less coherent than that for GRO 1655-40 which is 
presumably due to scattering in strong ambient wind/corona (see ST06 for details). 




